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Abstract 

We present a statistical mechanics description to study the ground state of quantum systems. In 
this approach, averages for the complete system are calculated over the non-interacting energy 
levels. Taking different interaction parameter, the particles of the system fall into non-interacting 
microstates, corresponding to different occupation probabilities for these energy levels. Using 
this novel thermodynamic interpretation we study the Hubbard model for the case of two elec- 
trons in two sites and for the half-filled band on a one-dimensional lattice. We show that the form 
of the entropy depends on the specific system considered. 
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1. Introduction 

Statistical Mechanics (SM) provides useful concepts to study systems with large number of 
particles. For example, based on standard SM, Edwards [1] proposed a thermodynamic descrip- 
tion of granular matter in which thermodynamic quantities are computed as fiat averages over 
configurations where the grains are static or jammed, leading to a definition of configurational 
temperature. A numerical diffusion-mobility experiment of a granular system has supported the 
Edwards' statistical ensemble idea [2]. Another example is the relation between entropy and the 
horizon area of a black hole [3|], which provides a new approach for studying black holes and 
quantum gravity theory. Furthermore, four laws of black hole mechanics can be demonstrated us- 
ing this thermodynamic description. The microscopic origin of the black hole entropy, originally 
calculated thermodynamically, has been explained from string theory. |0] 

Recently, Cejnar et al. [|7|] analyzed quantum phase transitions in finite systems [8] by defining 
an analog of the absolute temperature scale connected to the interaction parameter of the Hamil- 
tonian. And thus, they were capable of establishing a thermodynamic analogy for the quantum 
phase transition. However, they did not identify the correspondence with statistical mechanics 
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and consequently the new scenario opened by this microscopic analysis. This correspondence 
and these consequences are the goal of this paper. 

Here, we use tools developed in SM to study the ground-state of quantum systems. We ob- 
serve that, for certain classes of quantum systems, taking different intensities of the interaction 
between particles of the system corresponds to taking different occupation probabilities for non- 
interacting microstates energy levels. With this observation we can define an analog of the abso- 
lute temperature scale in such a manner that it is possible to make a thermodynamic interpretation 
for the interaction in the ground-state of quantum systems. The Hubbard Hamiltonian is a typical 
model in which this approach can be applied. Here, we analyze two exact solvable limits of the 
Hubbard model. 

This paper is organized as follow. The formalism is described in Sec. [2] The study of the two 
exact solvable problems based on the Hubbard model are presented in Sec. [3] Finally, we present 
the conclusions in Sec. [4] 

2. Formalism 

The scheme of our formalism can be applied to a broad class of Hamiltonians defined as 



where we assume that Ho is a one-particle Hamiltonian operator and the interaction term is given 
by the V operator and T is the dimensionless interaction parameter. Here, we must consider that 
T > and the operator V is positively defined. In this way we have established that the energy is 
a concave function of T 

A good example of this class of Hamiltonians is the one of the Hubbard model |5]. In this 
model, which is amongst the most important magnetic ones, the eigenstates of the Hamiltonian 
in the absence of interaction (T = 0) are just the non-interacting states |0,-), whose respective 
energy eigenvalues, £,(0) are defined through the relation Ho\<pi) - Ej{Q)\<p,). The eigenvalues 
Ei(T) of H for nonvanishing T are obtained from the equation H\xjJi) — Ei(T)\\jji). Moreover, the 
expectation value of an operator X on the ground-state |^o) is given by (^qIa^o). 

Now, we can provide a approach for obtaining expectation values of physical quantities on 
the ground-state in the base of non-interacting states. This simply means to find the expectation 
values of X on the ground-state in the representation. 

The ground-state can be expanded in terms of the non-interacting states \<pi) as 



where the coefficients aff) = (0,|t/<o)- We recall that the quantity \ai(T)\ 2 has a probabilistic 
interpretation. In other words, we can write pt(T) = |a,(r)| 2 e [0, 1] and Pi(J) = |a,(r)| 2 = 
1 . This establishes the connection to SM: the expectation value of Hq, 



can be interpreted as a usual average, which is computed over the set of non-interacting energy 
levels Ej(0), each one with probability pi(T). It is an analog of the mean energy (E) = 
We can easily verify that if T > and (V(T)) > is a monotonically decreasing function of T, 
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Figure 1: Similar to the usual canonical ensemble of the SM, taking different ground-state temperature T g (interac- 
tion parameter), the particles of the system fall into non-interacting microstates, corresponding to different occupation 
probabilities for these energy levels. 



then Eo(T\) < £0(^2) for T\ < Ti- In this case, in analogy to SM, for the non-interacting case 
T — 0, the system has the lowest energy Eq(Q) and pt(0) = (5,o- If T > 0, like a thermal energy, 
the interaction favors other energy levels of the non-interacting case. In this description, only 
the non-interacting microscopic states are used to compute the thermodynamic properties. This 
enables us to define an analog of the absolute temperature scale, called ground-state temperature, 
as T g = T jk, where k is a constant measured in Kelvins -1 . 

This description is illustrated in Fig. [T] Similar to the usual canonical ensemble of the SM, we 
can consider that taking different ground-state temperatures T g , i.e, different values of the inter- 
action parameter, the particles of the system fall into non-interacting microstates, corresponding 
to different occupation probabilities for these energy levels. 

In addition, an analogy to the standard thermodynamics is also reproduced by this description. 
We can introduce a so-called ground-state thermodynamics, defining the ground-state internal 
energy, ground-state free energy and ground-state entropy, respectively, as 

U(T g ) = (E(T g )) = £ Pi (T g )Em, (4) 

i 

F(T g ) = E Q (T g ) - r,<V(0)>, (5) 

S(T g ) = k((V(0))-(V(T g )}). (6) 

It can be easily seen that S(T g ) is a non-negative monotonically increasing function of T g . We 
can trivially verify that, using Eqs. ©-(O, the ground-state thermodynamics precisely satisfies 
the standard thermodynamics relation for the Helmholtz free energy 

F(T g ) = U(T g )-T g S(T g ). (7) 

Furthermore, we can derive the thermal response function, in correspondence to the heat capacity 

dS{T g ) d 2 F(T g ) 
C(T g ) = T g ^f = ~T g ^f • (8) 

3 



It is interesting to observe that the expression above can be calculated using the Hellmann- 
Feynman theorem which allows to find the ground-state expectation values of a general operator 
X by differentiating the ground state energy of a perturbed Hamiltonian Hq + AX with respect to 
A Bl 



3. Applications 

For illustrating the approach introduced in this letter, let us study two exact solvable problems 
based on the Hubbard model 051]. The Hamiltonian of the Hubbard model is defined as 

H = -t^ c] a c ja + Uj^ %%> (9) 
<y>« ' 

where c! , c la and = c ! c /Q . are respectively the creation, annihilation and number operators for 
an electron with spin a in an orbital localized at site i on a lattice of N sites; the (ij) denotes pairs 
i, j of nearest-neighbor sites on the lattice; U is the Coulombian repulsion that operates when 
the two electrons occupy the same site; and t is the electron transfer integral connecting states 
localized on nearest-neighbor sites. First and second terms of Eq. (0 correspond to, respectively, 
one-particle Hq and interaction terms of Eq. ([TJ. 

The problem of two electrons in two sites is the simplest example to our approach. By using 
direct calculus, it is easy to obtain the ground-state eigenvalue and eigenfunction, respectively, 
as 

Eo(U) = - 1 -{U - Vf/ 2 + (4f) 2 ), (10) 

and 

l«Ao> = a-\<(>-) + a + \<f> + ). (11) 
where \<p ± ) are eigenfunctions for the case U — 0, with a_ = -yl - a\ and 



a + = 2t/ y](2 ^U 2 + (At) 2 - U) ^U 2 + (4fp. (12) 

Thus, we define T g = U/kt, and using Eqs. (O-© into Eqs. (TTOt and (fTTT i. we find (from now 
k — 1 and t — 1 for simplicity) 

F(T g ) = -1 Jrj + 16, (13) 

S(T g ) = — (14) 
2^716 

U(T g ) = -J , (15) 



and 



87. 



C ^ - wuwr- (16) 



Figure|2]shows curves (full lines) of F(T g ), U(T g ), S(T g ) and C(T g ) versus the temperature T g for 
the expression above representing the case of two electrons on two sites for the Hubbard model. 
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Figure 3: Functional forms of the entropies S(p) associated with the Boltzmann-Gibbs (dashed line) and with Eq. U8t 
for 2 states (full line). 

As clearly seen in these figures, the behavior of these new variables is exactly as expected from 
the usual thermodynamics. 

Now, let us consider the functional dependence for the entropy given by Eq. ( fT4t in terms 
of the occupation probability of the non-interacting quantum states. Using the energetic con- 
straint (Eq. (H|), this dependence generates the concept of thermostat temperature, if we focus 
on the canonical ensemble of the SM formalism. It is easy to show from Eqs. (fTTb - (fT2l > that the 
occupation probabilities for the eigenfunctions \(f> ± ) of the non-interacting case are 

P± (T g ) = \ T 2 . (17) 

2 V r * 2+16 

We straightforwardly obtain the entropic form 

S{p)=yfp^pZ, (18) 

which is a concave function representing the geometric average of the quantum states probability, 
where certainty corresponds to S = 0. Here, we can see the difference between the standard SM 
and the ground-state SM. For the standard SM we always use the Boltzmann-Gibbs entropy 
S (p) = 2, pi In pi, while for the ground-state SM, this universality is broken, and the form of the 
entropy depends on the particular quantum system. On the other hand, this issue does not rule 
out the possibility that many different systems fall into some basic classes exhibiting qualitatively 
similar behavior. In Fig. [3] we show the functional forms of the entropies associated with the 
Boltzmann-Gibbs and with Eq. ( fT8l assuming 2 states. 

In what follows, we shall illustrate the above procedure by addressing the exact solution for 
the half-filled band of the Hubbard model for the one-dimensional case in the thermodynamic 
limit. This famous solution was obtained by Lieb and Wu in the sixties [12] using the Bethe 
anzatz. Since then, this result is considered one of the most important ones, owing to the lack 
of exact solution for the Hubbard Model. The ground-state as a function of the electron-electron 
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interactions U, for N sites in the limit N — > oo, is given by 



E (U) = -4N r J ^ J ^ dw , (19) 
Jo w[l +exp(wf//2)]' 

where Jq(w) and J\(w) are Bessel functions. It is, then, simple to obtain the quantities associated 
to the ground-state thermostatistics: 



F(Tg)/N = 4- 4 f ^ W,Jl ) W> 7,^ (20) 
4 Jo 



Jq(w)J\ (w)dw 
w[l + exp(vi'r i ,/2)] ' 



5( r)/iv=i-i r (21) 



4 2 Jo cosh W 4 ) : 

. J(){w)J\(w)f(w, T R )dw 
U(T g )/N = - *' , (22) 



-f 

Jo 



w[l + exp(wr s /2)] 2 ' 
where /(w, T g ) = [4 + (4 + wr g /2) exp(w^/2)] and 

T q r°° wJ Q (w)Ji (w) sinh(wr„/4)c/w 



f 

Jo 



C(T g )/N=-l * , , „. ,, • (23) 



cosh J (wr„/4) 



We show the behavior of F(T g ), U(T g ), S(T g ) and C(T g ) versus the temperature T g for the so- 
lution of the one-dimensional Hubbard model in Fig. [2] These curves correspond to the dotted 
lines and, as well noticed from the case of two electrons, they are also expected from the usual 
thermodynamics. 



4. Conclusions 

In summary, we introduce an approach to solve problems of quantum mechanics using con- 
cepts of statistical mechanics. We can consider that taking different ground-state temperatures 
T g , i.e, different values of the interaction parameter, the particles of the system fall into non- 
interacting microstates, corresponding to different occupation probabilities for these energy lev- 
els. 

We found that the functional form of the ground-state entropy depends on the particular quan- 
tum system. The break down of universality of the entropy is consistent with the concept of 
generalized entropies associated with a specific quantum Hamiltonian. 

Finally, the ideas presented here can eventually provide a mechanism for new approximation 
methods, such as the usage of the geometric average of the quantum states probability in the 
high dimensional limit for the Hubbard model. We can envisage in further works the study of 
the possibility that many different systems may fall into some basic classes of the ground-state 
entropy. 
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